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Abstract 

The scaling properties of the free energy, the diamagnetic moment, and the 
diamagnetic susceptibility above the phase transition from the ferromagnetic 
phase to the phase of coexistence of ferromagnetic order and superconductiv- 
ity in unconventional ferromagnetic superconductors with spin-triplet (p-wave) 
electron paring are considered. The crossover from weak to strong magnetic in- 
duction is described for both quasi-2D (thin films) and 3D (bulk) superconduc- 
tors. The singularities of diamagnetic moment and diamagnetic susceptibility are 
dumped for large variations of the pressure and, hence, such singularities could 
hardly be observed in experiments. The results are obtained within Gaussian 
approximation on the basis of general theory of ferromagnetic superconductors 
with p-wave electron pairing. 
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1 Introduction 



The discovery of coexistence of ferromagnetism and bulk superconductivity in Uranium- 
based intermetallic compounds, UGe 2 [H [21 [3], URhGe [I], UCoGe |6] has led to 
renewed interest in the interrelationship between ferromagnetism and superconductiv- 
ity. In these itinerant ferromagnets, the phase transition to superconductivity state 
occurs in the domain of stability of ferromagnetic phase, including sub-domains, where 
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a considerable spontaneous ferromagnetic moment M is present. This seems to be 
a general feature of ferromagnetic superconductors with spin-triplet (p-wave) electron 
pairing [TJEIIS] (see also reviews [TUIITT]). In such situation the thermodynamic proper- 
ties near the phase transition line may differ from those known for the superconducting- 
to-normal metal transition. 

The basic thermodynamic properties of these systems are contained in their T — P 
phase diagrams. According to the experiments [H [21 0, SI El E] carried on the above 
mentioned compounds and ZrZn2 [T2] , the T — P diagrams exhibit several basic fea- 
tures, which are shown in Fig. [1] (as claimed in Ref. [13J, owing to a special treatment of 
sample surfaces, only a surface superconductivity has been proven in ZrZn2). As seen 
from Fig. [TJ the phase transition line Tp(P), corresponding to the phase transition 
from normal (paramagnetic) state (N) to ferromagnetic phase (FM) is substantially 
above the lines referring to phase transition from FM phase to the phase of coexistence 
of ferromagnetism and superconductivity (in short, FS phase, or, FS). The exception 
is only very near to the critical pressure P c , where both ferromagnetism and supercon- 
ductivity vanish and phase transition lines are close to each other. This picture reflects 
the real situation in the above mentioned compounds, for example, UGe2, where the 
Tp(P) at ambient pressure P a is of order 53 K, whereas the maximal Tps does not 
exceed 1.23 K; for UGe 2 , P c ~ 1.6 GPa. 

Moreover, Fig. [1] shows that the line Tp^(P) of FM-FS phase transition may have 
two or more distinct shapes. Beginning from the maximal (critical) pressure P c , this 
line may extend to all pressures P < P c , including the ambient pressure P a ; see the 
almost straight line containing the point 3 in Fig. [TJ A second possible form of this 
line, revealed by experiments, for example, in UGe2, is shown in Fig. CD by the curve 
which begins at P ~ P c , passes through the point 2, and terminates at some pressure 
Pi > P a , where the superconductivity vanishes. These are two qualitatively different 
physical pictures: (a) when the superconductivity survives up to ambient pressure, and 
(b) when the superconducting states are possible only at relatively high pressure (for 
UGe 2 , P x ~ 1 GPa). 

Within the general phenomenological theory of spin-triplet ferromagnetic superconduc- 
tors [7], these different pictures are distinguished by simple mathematical conditions 
on the theory parameters 0, Ej. Therefore, there are both experimental and theoretical 
arguments to classify these superconductors in "type I" and "type II" p-wave ferro- 
magnetic superconductors, as proposed for the first time in Ref. [8] and illustrated in 
Fig. [TJ above the respective T FS (P) lines. The tricritical points 1, 2 and 3 (see Fig. [TJ), 
at which the order of the phase transitions changes from second order (solid lines) 
to first order (dashed lines) are a quite reliable experimental fact and have already 
been theoretically explained [TJ. Firstly, the interaction between the superconducting 
and magnetic subsystems naturally generates a first order phase transition along the 
high-pressure part (P < P c ) of the FM-FS phase transition. Secondly, under certain 
circumstances, an additional M 6 term in the free energy may describe the experimen- 
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Figure 1: An illustration of T — P phase diagram of p-wave ferromagnetic supercon- 
ductors (details are omitted): N - normal phase, FM - ferromagnetic phase, FS - phase 
of coexistence of ferromagnetic order and superconductivity, T F (P) and T F s(P) are 
the respective phase transition lines: solid line corresponds to second order phase tran- 
sition, dashed lines correspond to first order phase transition; 1 and 2 are tricritical 
points; P c is the critical pressure, and the circle C surrounds a relatively small domain 
of high pressure and low temperature, where the phase diagram may have several forms 
depending on the particular substance. The line of the FM-FS phase transition may 
extend up to ambient pressure (type I ferromagnetic superconductors), or, may termi- 
nate at T = at some high pressure P = P\ (type II ferromagnetic superconductors, 
as indicated in the figure). 

tally observed first order phase transition along the high-pressure part of the N-FM 
line 0. 

In Fig. (U the circle C denotes a narrow domain around P c at relatively low temperatures 
(T < 300 niK), where the experimental data are quite few and we may not reliably 
conclude about the shape of the phase transition lines in this T — P domain. It could 
be assumed, as in the most part of the experimental papers, that (T = 0, P = P c ) is 
the zero temperature point at which both lines Tp{P) and Tps{P) terminate. A second 
possibility is that these lines may join in a single (N-FS) phase transition line at some 
point (T > 0, P' c < P c ) above the absolute zero. In this second variant, a direct N-FS 
phase transition occurs, although this option exists along a very small piece of N-FS 
phase transition line: from point (0, P c ) to point (T > 0,P C ' < P c ). A third variant is 
related with the possible splitting of point (0, P c ), so that the N-FM line terminates at 
(0, P c ), whereas the FM-FS line terminates at another zero temperature point (0, Po c )] 
P 0c < P c . In this case, the p-wave ferromagnetic superconductor has three points of 
quantum (zero temperature) phase transitions [9]. 

These and other possible shapes of T — P phase diagrams are described within the 
framework of the general theory of Ginzburg-Landau (GL) type [7j in an entire con- 
formity with the experimental data [HI E]; for reviews, see Refs. [TOT ITT] . The same 
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theory has been confirmed by a microscopic derivation based on a microscopic Hamil- 
tonian including a spin-generalized BCS term and an additional Heisenberg exchange 
term [14] . 

Although the theory predicts correctly the shape of FM-FS phase transition line, 
the possible types of phase transitions, finite temperature and quantum multi-critical 
points [9j [15], some important features of p-wave ferromagnetic superconductors, in 
particular, the properties of the FM-FS phase transition, need a further investigation. 
The FM-FS phase transition in p-wave superconductors is remarkable for the circum- 
stance that there the superconductivity appears in an environment of a strong ferro- 
magnetic moment (magnetization density) M. This certainly leads to a modification 
of the usual phase transition to superconducting state in conventional non-magnetic 
superconductors. 

In this paper we study the effect of the spontaneous ferromagnetic moment M deeply 
below the ferromagnetic transition on the fluctuation properties in a close vicinity of the 
FM-FS phase transition above the FM-FS line. Remember, that near the usual critical 
point of standard (conventional) superconductors, the fluctuating superconduction field 
ip(x) creates an overall diamagnetic moment VWldia = M. = = 1,2,3} which, 

except for particular circumstances, depends on the external magnetic field H. The 
numerous studies of the fluctuation diamagnetism in conventional superconductors, 
including various sample geometries [three-dimensional (3D) (bulk), quasi-2D (in short, 
q2D; thin films), ID (wire), OD (small drop)] and layered structures, are summarized 
in Ref . [IE] ; D denotes the effective dimension of the superconductor places in a three 
dimensional space (d = 3); note, that the effective dimension of the superconductor D 
is often different from the space dimensionality d. 

Here we consider for the first time the diamagnetic properties of 3D and q2D p-wave 
ferromagnetic superconductors — diamagnetic moment Ai and diamagnetic suscepti- 
bility Xdia- Some of our preliminary results for 3D (bulk) superconductors have been 
recently published [17]. We demonstrate that the diamagnetic susceptibility Xdia and 
diamagnetic moment M. of p-wave ferromagnetic superconductors are dumped along 
the prevailing part of FM-FS phase transition line, in particular, at the second order 
phase transition line. So, our theoretical results predict that the singularities at crit- 
ical points, typical for usual superconductors, do not exist at the second order phase 
transition line Tps{P) of p-wave ferromagnetic superconductors. 

In Sec. II we present the GL fluctuation Hamiltonian of p-wave ferromagnetic supercon- 
ductors in Gaussian approximation for the fluctuation field ^(x) of the superconducting 
order parameter ip. We consider temperature T > Tps(P), where the statistical av- 
erage (VK X )) of if) = (^( x )) + Sif)(x) is equal to zero, namely, the field i/j{x) is a net 
fluctuation: ip = 6ip(x). The Gaussian approximation is the usual tool for study of 
basic properties of fluctuation diamagnetism [IB] , and here we follow this approach as 
well as the notations in Refs. [TJ El [TEJ US]- We shall essentially use formulations and 
results of the general phenomenological theory of p-wave ferromagnetic superconduc- 
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tors [7J El El QUI EH] , as well as modern concepts of crossover and critical phenomena [TS] . 
Our theoretical approach essentially generalizes preceding theoretical treatments (see, 
e.g., Ref. [16J and references therein) and might be used as an advanced calculational 
scheme in further investigations of critical crossovers in complex superconductors. 

In Sec. Ill we derive a general expression for the fluctuation contribution to the equi- 
librium free energy of system above the curve Tps(P) for arbitrary values of magnetic 
induction B = H + 47rM. We show that apart of special role of magnetic induction B, 
this expression is very similar to that for conventional nonmagnetic superconductors. 

In Sec. IV we calculate the diamagnetic moment M.(T, H) and the diamagnetic suscep- 
tibility Xdia{T, H). Further, we investigate the "weak-i? - strong- B n crossover in the 
behavior of these important quantities. For this aim, both 3D and q2D sample geome- 
tries are considered and the results are compared with known studies of conventional 
(s-wave) non- ferromagnetic superconductors. We predict a dumping of the singulari- 
ties of the diamagnetic moment and the diamagnetic susceptibility along the FM-FS 
phase transition line. Moreover, we deduce a universality in the behavior of these 
quantities in p-wave ferromagnetic superconductors and usual (s-wave) non-magnetic 
superconductors. We will also briefly focus on related experimental problems. In Sec. 
IV. C and Sec. V we summarize our main results and discuss their applicability to real 
systems. 



2 Fluctuation Hamiltonian 



The fluctuation Hamiltonian of p-wave ferromagnetic superconductor can be given in 
the form [7j 



2 

+ a s \i)\ 2 

+ yH 4 +*7oM. (V> x r) + 5oMV| 2 ] , (1) 

where 2m and 2e = 2|e| are the effective mass and the charge of the electron Cooper 
pairs, respectively, the superconducting order parameter ip = {ipj(x);j = 1,2,3} is 
a three-component complex field, the magnetization M(x) = {Mj(x);j = 1,2,3}, de- 
scribing the ferromagnetic order in the FM phase is a three component real field, i.e., 
the components M?(x) are real fields, a s = a s (T — T s ) is represented by the generic crit- 
ical temperature T s of hypothetical pure superconducting state (|^| > 0, M = 0) and 
the positive material parameter a s ; as usual, b s > 0. Besides, 70 ~ J > 0, where J > 
is an effective ferromagnetic exchange constant, and 60 are parameters describing the 
interaction betwee superconducting and magnetic electron subsystems. As usual, the 
vector potential A = {A,(x)} obeys the Coulomb gauge V.A = and is related to the 
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magnetic induction B. In a space of dimensionality d, in short, dD dimensional space, 
the relation between A and B is represented in different mathematical forms. Here we 
work in a three dimensional space (d=3), where the geometry of superconductor body 
is chosen of two types: q2D (thin film) or 3D (bulk); hence, we may use the relation 
B = V x A. We neglect the gradient anisotropy [20] as it has a small effect compared 
to the exchange interactions between the normal and superconducting electrons; this 
point is discussed in the Section below. 

In mean field (MF) approximation [T8| [19] . the magnetization M = |M| in the pure 
ferromagnetic phase (FM) (ip = 0, M ^ 0) at zero external magnetic field (H = 0) 
is given by Mmf = (~ a f /bf) 1 ^ 2 as a global minimum of the magnetic fluctuation 
Hamiltonian 



n{M) = J d 3 x 



. 3=1 



(2) 



where a/ = a/(T — T F ), T F is the generic critical temperature of the (pure) ferromag- 
netic state, and «/ and bf are positive material parameters [7]. Note, that throughout 
this paper M denotes the density of the magnetization, whereas the total magnetization 
of the ferromagnetic phase is given by / d 3 iM(x). 

In the present investigation we ignore the fluctuations of the field M(x). This is 
justified by the fact, that T F (P) 3> T FS for almost all variations of the pressure P, 
except for a very narrow domain near the critical pressure P c , where the mean-field 
value M MF = (— a//6/) 1//2 of the magnetization tends to zero (T FS ~ T F ). Remind, in 
this domain both N-FM and FM-FS phase transitions are of first order, and therefore 
strong fluctuations of M could not be expected. 

The last two terms in Eq. (CD) have a key role in the description of thermodynamics and 
phase diagram of this type of ferromagnetic superconductors. These terms describe 
the interaction between the normal (non-superconducting) electron fraction and the 
electron Cooper pairs. While the 5o-term has the supporting role of ensuring a stability 
of phases at relatively large negative values of the parameter df, the 70-term has a key 
role in the phenomenon of coexistence of superconductivity and ferromagnetism. The 
presence of this term ensures the theoretical description of real situation in p-wave 
ferromagnetic superconductors, in particular, a reliable description of the mentioned 
coexistence of phases. This is the term which triggers the superconductivity [?]; for 
a more detailed discussion, see Refs. [7J [10j [H]. This point is not trivial. Although 
such "trigger" terms, where one of ordering fields is present by its first power and the 
second order parameter interacts by its second power, are well known, for example, in 
the theory of certain improper ferroelectrics [21], here the symmetry of the 70-term is 
very particular and leads to critical phenomena of essentially new universality class |1 5]. 



Note, that according to the symmetry analysis in Ref. [22J, a 70-term follows from the 
gradient anisotropy, which is typical for unconventional superconductors [20]. How- 
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ever, if the gradient anisotropy is the only source of creating such important term, the 
thermodynamics of these systems could not be described in compliance with the exper- 
imental data. This is so, because the exchange energy that creates the ferromagnetic 
order is of much bigger magnitude than the potential corresponding to the occurrence 
of FS. This argument is readily justified by the inequalities T F 3> T FS 3> T s ~ 0, 
which follow from the experimental data. Therefore, the major contribution to the 70- 
term and to the effective interaction parameter 70, comes from the interaction between 
the Cooper pair fraction and the normal (non-superconducting) electrons of conduc- 
tion electron bands of compounds rather than from the exchange interaction between 
Cooper pairs in the superconducting sub-band only. So, contrary to the consideration 
in Ref. [22], where only the superconducting sub-band takes part in the description 
and 70-term is a net product of the gradient anisotropy of p-wave Cooper pairs, within 
the present general quasi-phenomenological approach, the interaction parameter 70 in- 
cludes both the inter-sub-band exchange interaction between the magnetic moments of 
Cooper pairs and normal electrons, and the exchange interactions within the sub-band 
of Cooper pairs, namely, between the Cooper pairs themselves. Obviously, the second 
type of interaction is much weaker, and can be safely ignored in many calculations. 
This point of view is supported by recent microscopic theories [MJ 123] . 

As known, for the Uranium-based compounds, all these remarkable phenomena are 
created by the 5/-electrons, whereas for ZrZn 2 the same phenomena might be ascribed 
to the behavior of the 4<i-electrons. Furthermore, within our approach one may extend 
the consideration beyond the itinerant ferromagnetism and take into account exchange 
effects on the conduction band electrons produced by localized spins, attached to atoms 
at the vertices of crystal lattice. 



3 Free energy 

Here our task is to calculate the fluctuation part of free energy 

3 

F = -/T 1 In / J] II ^i( x ) WM iW ex P ( 3 ) 
J iJ=lxeV 

in the superconductor volume V = L x L y L z , above the FM-FS phase transition line 
Tfs{P) = T C (P); henceforth we shall use the notation T c for the temperature Tpg- In 
Eq. ([3]), the functional integral is over all independent degrees of freedom: the complex 
components of the superconducting fluctuations, i/j(x) and the fluctuations <5Mj(x) of 
the magnetization vector components Mf, = k^T . The functional integral is taken 
over both real [9fJ?/>(x)] and imaginary [Q^(x)] parts of the complex field ^(x), i.e., 
2>0(x) = c$Fh/>(x)cZ9^0(x). Note that for temperatures near T C (P) we can always set 
(3 ?s (3 C = 1/ksTc (see, e.g., Ref. [18J). Mostly in this paper we shall consider 3D (bulk) 
superconductors and (q2D) thin superconducting films in a transverse magnetic field, 
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so we shall use the 3D notations: x = (x,y,z), and the labels x,y,z for quantities 
defined along the respective Cartesian axes in 3D space (d — 3). 

We are interested in the magnetic thermodynamics in close vicinity (T ~ T c ) above 
the temperature T C (P) of FM-FS phase transition. For Tp 3> T c , the fluctuations 8Mi 
in this temperature range are very weak and can be neglected [M< « M^ (MF) ]. Then 
one may substitute Mj(x) in the exponent of Eq. (E]) by Mj Mi?) and V8Mj( x) in the 
functional integration by T , 5[5M 7 (x)]; here 5 [2] denotes 5-function. This procedure 
totally eliminates the fluctuations <5Mj(x) from our consideration; hereafter we shall 
omit the label "MF" [M( MF ) = M]. For T > T FS the statistical averages (ipj) = 
■ipj — 8ipj are equal to zero and we have a purely fluctuation field ip(x) = 6if)(x). Thus, 
within the present consideration, the only integration variables in functional integral ([3]) 
are the fluctuations 8ipj{x) = Neglecting the Ginzburg critical region [HI H~9] . 

which is very small in all low temperature superconductors and, hence, unobservable 
in experiments, we may ignore the fourth order fluctuation term \ipj\ 4 in Eq. (JTJ, and 
apply Gaussian approximation to the fluctuation modes VK X )- 

For convenience, we choose the vectors M and H along the z-axis: M = (0,0, M), 
and H = (0,0, if). This assumption does not essentially restrict the generality of 
our consideration. If the superconductors is magnetically isotropic, the magnetization 
vector M will follow the direction of external magnetic field H. If magnetic anisotropy 
is present, for example, an easy axis of magnetization, as in UGe2, our assumption will 
be satisfied by choosing an external field H parallel to this easy axis of magnetization. 
Then the term M.(-0 x ifj*) takes the simple form M(ipiip2 — c.c). 

Under the supposition of uniform magnetic induction B = (0, 0, B), we take the gauge 
of the vector potential A as A = (—By, 0, 0) and expand the fields %(x) in series [19] 

^•( x ) = l^ c j(g)^-(g, x ) ( 4 ) 
1 

in terms of the complete set of eigenfunctions 

Pj(g,x)= ) , 1/2 e t(W Xn(?/) (5) 
of the operator [ihV + (2e/c)A] 2 /4m, corresponding to the eigenvalues 

E(q)= ( n +0^+|^> ( 6 ) 

specified by the magnetic frequency ujb = (eB/mc) and vector q = (n,k x ,k z ), where 
n = 0,1,.. .,00, is the quantum number corresponding to the Landau levels, and k x 
and k z are components of the wave vector k = (k x , k y , k z ). In Eq. (jSJ), the function 
Xn(y) is related to the Hermite polynomials H n (y) by 
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Xn (y) = A n e >«b H n , (7) 

V a B J 

where A' 1 = (a B 2 n n\ v ^) 1 / 2 [24J, y = and a B = (hc/2eB) 1 ' 2 . 

In terms of the Cj(q) -functions, the ^> 2 -part of the fluctuation Hamiltonian (CQ) is given 
by 

ft = J2 E(q)cj(q)c*(q) + i lo M [ Cl (g) C *(g) - c.c] , (8) 



3,1 



where 



E(q) =E(q) + a s + S M 2 . (9) 
Applying the unitary transformation, 

Cl (g) = -J=H> + (g) + <Mg)] (10a) 

c 2 (g) = i=[0 + (g) + <Mg)] (10b) 

renders the fluctuation Hamiltonian (JTJ) as a sum of squares of field components c 3 (g), 
and <f)±(q)- 

In the continuum limit, L y oo, the integral 

/Ly/2 
dyxn(y)xn'(y) (n) 
-L y /2 

is simply equal to the Kronecker symbol 5 nn > and this is a key point in the further simple 
representation of the Hamiltonian ([1]). Substituting the function (J7]) in the integral 
( ITT]) , and having in mind the properties of the Hermite polynomials H n (z) [23], we 
see that the integral (fTTj) will be equal to only if the limits of integration can be 
expanded to ±oo. In fact, the integration in Eq. ( flT]) can be performed with respect 
to the variable y = (y — yo)/ a £ and integral limits (— yo ± L y /2)/ctB- The limits of 
integration with respect to this variable can be approximately equalized to ±oo, if only 
L y /2as — > oo; for finite samples (L y ^> 2a^). Another inevitable condition, namely, 
—Ly/2 < y < L y /2, follows from the requirement that the coordinate y = a 2 B k x must 
belong to the sample volume. This condition implies 

Ly < k < h-- (12) 



2a| x 2a| 



As we show below, this condition fixes the number of states M = L x L y /2na 2 B for all 
possible values of quantum number k x at any given n and k z . 
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Having in mind these features of theory and considering sufficiently large L y we can 
justify the solution I y = 8 nn i of the integral ( ITTj) and achieve a very useful form of the 
Hamiltonian, namely, 



n=J2[E4q)\Mq)\ 2 + £ + (?)l<M?)l I 

n,q 

+£ 3 (g)|c 3 (g)| 2 ], 



with 



E±(q) = E(q) + a±{M), 
where E(n,k z ) is given by Eq. flH]), 



(13) 
(14) 



a±(M) = a ±7 M 



(15) 



is represented by a = a s + S M 2 , and E 3 = E is given by Eq. (jUJ). 

Now the free energy (j3J) can be written as a functional integral over all independent 
field amplitudes: ±0(g), and 03(g). Using the short notations <^ Q (g) with a = (+, — , 3) 
of the Fourier amplitudes <f>+(q), 4>-{q) and 03(g), respectively, and adopting the same 
label a to denote E + , EL and £"3 by E a , we obtain the free energy F in the form 



The direct calculation of the Gaussian integrals in Eq. ffl6|) yields 



(16) 



F = -k B T In J] 
= /. /; 71n J"} 



EJn, k z ) 



TlksT 

E a (n,k z ) 



(17) 



where we have used the condition ffl2|) and the continuum limit for the /c^-product, 



exp 



L x J dk x /2n 



namely, 



L y' 2a2 B dk x L V L X 



—L y /2a- 



2 2/T 2ira 



2 ■ 



(19) 



In the second equality ( TT71) we point out the result of the summation over k x with the 
help of the rule (jl8jl . The obtained expression (fl71) follows from the fact that the mode 
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energies E a (n,k z ) do not depend on the quantum number k x . Thus the number of 
states M at fixed quantum numbers n and k z and the relation F ~ J\f are naturally 
deduced from the calculation. 

Further, we have to pay attention to the fact that the field theories of GL type are lim- 
ited to length scales k = |k| < A ~ 7r/£ , where £o is the zero temperature correlation 
length of the field of interest [TH [19]. In the present case, we must use the quantity £ 
corresponding to the field ip(q) which fluctuates in the vicinity of the phase transition 
line T C (P). The standard expression £q s = h/ 1 (Ama s T s ) 1 / 2 \ corresponding to the generic 
critical temperature T s [191 IB] cannot be applied to our problem. So, we define the 
upper cutoff A ~ 7r/£ by the zero-temperature correlation (coherence) length £ but 
the latter will be specified at a next stage of our consideration. Here we will men- 
tion that £ is the scaling amplitude of the correlation length £ of the superconducting 



fluctuations at the FM-FS critical line T C (P): f (t) = Wl^| 1/2 , where t = {T - T c )/T c ; 



As the small wave numbers k have the main contribution to the values of the integrals 
in the free energy and its derivatives, we shall use the finite cutoff A only when the 
respective integral has an "ultraviolet" divergency; for example, see Eq. (fl9l) . In all 
other cases, the relatively large values of k do not produce essential quantitative con- 
tributions to the integrals and for this reason, we may extend the cutoff A to infinity. 
Moreover, owing to the same type of limitation of the GL theory — the long wave- 
length approximation £ k < tt, we should take in mind that only quantum numbers 
n corresponding to energies E a < S 2 /4m^ are to be taken into account. Thus the 
quantum number n has a cutoff as well, and the latter is given by ti c ujb — ft/4m^, 
namely, n c = [H/ 'Aitl^qUb] {[z] denotes the integer part of number z). This energy cutoff 
could be neglected in cases when this does not produce divergencies of the respective 
physical quantities. 

For our further aims we shall write Eq. (fTT]) in a more convenient form: 



The relevant part of the free energy, which contains singularities at the critical temper- 
ature T c is given by the term containing In EL. All other terms are quite smooth near 
the line T C (P) and do not produce singularities of the physical quantities. This impor- 
tant circumstance follows directly from the fact, that namely the parameter a_(M) is 
relatively small in magnitude and changes sign at T c — the FM-FS phase transition 
temperature, corresponding to the phase transition from FM phase to the phase of 
coexistence (FS) of ferromagnetism and a homogeneous (Meissner) superconducting 
state at zero external magnetic field (H = 0). Therefore, the critical fluctuations are 
described by the field 4>-{q). The other fields, 4>+(q) and 03(g), do not produce criti- 
cal phenomena (singularities) because the parameters a + and do not pass through 



t\ < 1. 




(20) 
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the null at T = T c . The value of magnetization M is relatively large along the most 
part of the line T C (P) in this type of ferromagnetic superconductors (Tp 3> T c ), and 
therefore the parameters a + and a$ are quite different from a_ except for a narrow do- 
main around the critical pressure P c . Note, that in a non-magnetic (standard) p-wave 
superconductor [18J, where M = 0, all three modes <p a {q) are critical 

in a close vicinity of the critical point T s . For such superconductor, the free energy 
( 12T]) will differ with a factor 3 from the standard result for a conventional (s-wave) 
superconductor with a scalar order parameter [19] . 



In our further analysis we shall ignore the nonsingular part of the free energy and keep 
only the contributions from the critical mode (j)-(q). Thus we have to analyze the 
behavior of function 

F = - S e -^I1 f y In - ^ bT ,^ (21) 
where E(q) is given by Eq. ([6]). 

Now we have to define the parameters in Eq. (|2ip . For M = (|a/|/&/) 1/ ' 2 , the parameters 
a±(M) are given by [HI E] 

> 1/2 

fl / JL „, a f 



a ± (T) = a s (T-T s ) + 5 ^-± l0 \^J • (22) 

We are interested mainly on the parameter a_ (T) which is related with the equilibrium 
phase transition from FM to FS. Defining T c = Tps from the equation a_(T c ) = 0, we 
obtain 



a_(T)^a c (T-T c ), (23) 

where 

r I/ 2 

a c = a, - J r^ + 1 ( 24 ) 

h 2[b f {T F -T c )\ 1 ' 2 

and T c is given as a solution of the equation 



T C (M C ) = T s - ^M 2 C + ^M c , (25) 

1 /2 

where M c = M(T C ) = [af[Tp — T c )/bf\ > 0. In the same way one obtains that the 
parameters a and a + remain positive at T c : a Q (T c ) = 7oM c , and a + (T c ) = 2ao(T c ), 
which is a demonstration that the modes <f>+(q) and 03(g) are not critical and could not 
have essential contributions to the thermodynamics in the vicinity of phase transition 
line T C (P). 

Note that the parameter a c , given by Eq. ( 12"%|) . is positive for requirements of stability 
of the ordered phases [TD] • The solution of Eq. (125]) with respect to T c yields the curve 
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Tps(P) = T C (P), shown in Fig. [TJ The dependence of T c on the pressure P comes from 
the P-dependence of material parameters (a s , a/, Tp, . . . ) in Eq. ( )25|) . In Refs. [HI E] a 
simple P-dependence of these parameters has bee suggested: all material parameters 
except Tp are P-independent, and the form of the function Tp(P) is assumed of the 
simple form Tp{P) m Tp(0)(l — P/P c ). Although this is a simple approximation of the 
pressure effect in these systems, it gives a remarkable agreement between theory and 
experimental data for the T — P phase diagram [EJ |9]. In the framework of the same 
approximation, according to Eq. ( 1251) . the P-dependence is contained in the quantity 
M c and, hence, we may often consider T c as a function of M c = M(T C ) - the value of 
the magnetization on the FM-FS phase transition line: T c = T C (M C ). 

The upper cutoff for the wave number k z is given by A ~ tt/C,o, where the zero tem- 
perature correlation length £ — (h 2 / 4:ma c T c y/ 2 is expressed by a c and T c , given by 
Eqs. (I21I) and (|25|) . respectively. The upper quantum number n c , defined by the equal- 
ity n c hjjj c ~ h 2 /AmC,Q, can be represented by n c = [1/26], where b = [efiB /2mca c T c ) 
is a non-negative quantity. Having in mind the supplementary condition that E(n, k z ) 
from Eq. ([6]) with n = k z = should also obey the condition P(0,0) < h 2 /Am^, 
we find that b G [0, 1]. For type II superconductors, the parameter b has the useful 
representation b = P/P c2 (0) by the upper critical induction P c2 (T) = P C 2(0)|t| at zero 
temperature, P C 2(0) = eh/2mca c T c , i.e., for \t\ = 1; henceforth we shall denote B c2 (0) 
byP - 

Using these remarks, we can represent the free energy in the form 



F = pSB f(b,t) (26) 

where p = (eksT /Tihc), and 

A 

f(t,b)= S(k z ,t,b) (27) 

k z =-A 

is given by the sum 

n=0 u 2 

The function f(t, b) describes the shape of the free energy P(T, P), whereas the func- 
tions m(t,b) = —df/db and m'(t,b) = dm/db represent the variations of the dia- 
magnetic moment M. = —dF(T,H)/dH and the diamagnetic susceptibility Xdia = 
—d 2 F(T,H)/dH 2 , respectively. For the choice H = (0,0, H), the diamagnetic vec- 
tor Ai has only one component: (0,0, M). As OH = <9P, we can use the formulae 
M = —dF(T, B)/dB and X dia = dM(T, B)/dB. 

The wave number k z lies in the reduced Brillouin zone: — 7r/£o < k z = 2tcI/L z < 7r/£o; 
/ = 0, ±1,... ± [L z /£o\. For q2D systems, where L z < £ , the only possible value 
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of k z is zero and, hence, for q2D superconductors the function f(t, b) coincides with 
S(k z ,t,b) = S(0,t,b). For 3D systems, we shall use the continuum limit for f(t,b), 
given by 

/ A rlk 
^S(k„t,b). (29) 

In Eqs. (12 ip and ( |28l) . the logarithmic divergence at maximal temperature corresponds 
to n = 0, k z = 0, and e = (b + 1) — 0. The parameter t indicates the vicinity to T c 
along the T-axis and the parameter b = B/Bq shows the strength of the induction B 
and the distance to the phase transition point (T c , B = M) along the if-axis of the 
(T, P, H) phase diagram. These two parameters, t and b, are suitable for investigations 
of the system properties for t > 0. It is easy to show that for t < 0, where the upper 
critical induction £? c2 (T) = B (—t) > of type II superconductors is defined, the 
parameter e = (t + b) can be represented in the suitable form e = [B — _B c2 (T)]/£> , or, 
alternatively, in the form e = [T — T C 2(B)]/T C , where T C 2(B) = T c (l — b) is the higher 
critical temperature of type II superconductors [19]. Thus, for t < 0, the parameter 
e shows the distance from the phase transition line described by the upper critical 
induction B c2 (T) = B \t\, or, alternatively, by the higher critical temperature T c2 (B). 
The parameter e is appropriate for investigations at temperatures T < T C (P), i.e., 
t < 0. In this paper our consideration is restricted to temperatures T > T C (P) and for 
this reason we shall use the original parameters variables t and b, as given in Eqs. ( 1261) 
- (EH}. 

Performing the summation in Eq. ( )28l) and keeping only terms which depend on b ~ B, 
we obtain 



S(k z , t, b) = -b In ^ J - b ln(l + 1) 
+ ^ + 61n [l + b + e(k z )] 



r 


"1 + b + e(k z )' 


26 


r 


~b + e(k z )~ 




26 



(30) 



where e(k z ) = (t + £^k 2 z ) and T(z) is the gamma function. The sum S(k z ,t,b), given 
by Eq. (13"0"|) and the shape function f(t, b), given by Eq. (|2T|) . contain redundant terms. 
Remind that we have neglected the contributions from the factors E + and E 3 in Eq. ( 120]) 
although most of them depend on the parameter b, namely, on the induction B and, 
hence, these terms may have a finite contribution to the diamagnetic moment Ai. The 
mentioned terms have been however ignored for the fact that they do not produce 
singularities in the free energy derivatives. Thus, within the approximations already 
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made, we cannot evaluate correctly the magnitude of finite contributions to important 
quantities as M. and Xdia at the phase transition point. We may just demonstrate that 
such contributions exist. To be in a consistency with our preceding consideration, we 
should neglect such terms in Eqs. fl26|) . f )27|) and f )30|) . too. 

Up to now we keep all 6-dependent terms contained in the general formula ff2TT) for 
the free energy, including terms which obviously does not lead to any singularities, 
for example, the first and second terms on the r.h.s. of Eq. fl30|) . At this stage we 
make the stipulation to keep these terms in our further consideration with the remark 
that contributions to Ai and Xdia which are finite at the phase transition point will be 
neglected, provided divergent term is present. When no divergency occurs in some of 
this quantities, we shall keep the finite term only to indicate the lack of divergencies 
and to show that the respective quantity remains finite at the phase transition point. 



4 Crossover from weak to strong magnetic induc- 
tion 



The free energy, the diamagnetic moment A4 and the diamagnetic susceptibility Xdia 
can be investigated analytically in the limiting cases of strong and weak magnetic 
induction B. Above the critical temperature T c , when t > 0, the weak-5 limit is 
defined by the condition 

e(k g ) = t + ilk\ > b, (31a) 
whereas the strong- limit is given by the opposite condition 



e{k z ) < b. (31b) 

These conditions are considered in a close vicinity (t <C 1) of the phase transition line 
T C (P). The condition (13 lap is satisfied for any k z , provided t S> b. The condition 
( 131pp . however, can be satisfied only for t <C 6 and (£o^z) 2 sufficiently small. For large 
k z ~ tt/Co this condition does not hold for any b G [0,1]. However, the sum (12 7p 
is practically taken for k z ^ only for the 3D geometry, and in this case, the main 
contribution to the sum is given by the relatively small wave numbers {^k z ^ 7l )- 
Having in mind this argument, we can use the condition ( 131b|) without any restriction 
to small values of ^k z because the final result for the free energy will not essentially 
depend on the contribution of relatively large wave numbers. In particular, this is true 
in the continuum limit (I29P for the sum (12 7p and (t + b) 1, namely, in the close 
vicinity of the phase transition line T C (P). 
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4.1 Weak-5 limit 



Applying the condition (13 lap to the sum ( 1301) we obtain the result 



S{k z ,t,b) =b In 



e(k) 



6 2 
+ 12 



Vfc B /a c )(l + 1) 
11 1 



(32) 



In Eq. ( 1321) . all 6- independent terms have been omitted as irrelevant to our consideration 
and small terms of type 0(b 3 ) have been neglected. From Eq. (1321) is readily seen 
that only the term of type b 2 /e(k z ) will produce a singularity of the thermodynamic 
functions as this term tends to infinity for k z ~ and t ~ 0. Therefore, in this case, 
we can neglect the other terms in Eq. (132]) and write the free energy in the form 



P SB< 

125 n 



e{k g )' 



(33) 



This result has been obtained in Ref. [T7] in different notations. 



4.1.1 3D superconductors 

For 3D superconductors and tc 1, the sum ( 1271) over k z can be substituted by the 
integral ( 1291) and the cutoff A can be extended to infinity. Then the free energy becomes 

F - pVB2 CU) 

The diamagnetic moment Ai = M^ ia = —dF(T, B)/dB takes the form 

*»< r ' Jr > = -l2&- (35) 



In contrast to the usual case of non-magnetic superconductors [16] , here the diamag- 
netic moment does not vanish at H = but rather remains proportional to the magne- 
tization M c at the FM-FS phase transition line. In fact, as we work at a close vicinity 
of FM-FS phase transition line T C (P), the magnetic induction in Eq. (135|) should be 
approximated by B w B c = H + 47rM c , where M c = M(T C ). Within the weak- 1? limit, 
this result is valid for relatively small values of M c : M c B . 

The diamagnetic susceptibility Xdia^ = 9M-3D/9B is given by 

X ^ )(T) = - 1^ ' (36) 
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In these general notations this is the well known result for the fluctuation diamagnetic 
susceptibility above the critical point T c of conventional superconductors piJ[T9]. For 
p-wave ferromagnetic superconductors we have to take into account that T c = Tp$ as 
given by Eq. (125]) and the material parameter a c , which enters in the zero-temperature 
correlation length £ > is given by Eq. (I24|) . Thus one reveals the result for Xdia{T), 
obtained in Ref. |17j . 



4.1.2 Quasi-2D superconductors 

For thin films, where L z < £o = 7r/A, only the wave number k z = satisfies the con- 
dition — 7r/£o < k z = 2-kI/Lz < 7r/£ ; I = 0, ±1, . . . , [kL z /£q]. For such q2D geometry, 
f{t,b) = S(0,b,t). In the weak-£> limit (I31al) . S(0,b, t) is obtained by setting k z = 
in Eq. (1321) . Once again we may keep only the leading singular term b 2 /12t. Thus we 
obtain the free energy in the form 

F - pSB2 (Vf\ 

Fw ~ vim- (37) 

In Eq. f l37|) and below we use the label "2D" to denote quantities corresponding to q2D 
systems. Now one easily finds that 



« H > = -5i| (38) 

x2f»(T) = (39) 

Having in mind the relations f l2"Ej) and (|2"T|) as well as d/db = B d/dB, for q2D systems 
we obtain Xdia = (pS/B Q )x(t,b), where is the susceptibility shape function. 

The latter is defined by x = -d 2 f(0,t,b) /db 2 with f(t, b) = S(0,t,b), where S{0,t,b) 
denotes S(0,t, b) with {jikB/oi c ) = 1 for suitable choice of units, as given by Eq. ( 130|) 
for k z = 0. The function x(t,b) is depicted in Fig. [2] for t = (0.01, 0.1) and t = 
(0.01, ...,0.1). As shown in Fig. [2} the shape function x(t, b) exhibits a sharp decrease 
even at values of (t, b) ~ (0.01, 0.05. The minimal value x ~ — 10 for t ~ 6 ~ 0.01 is a 
precursor of divergency, given by Eq. As we see from Fig. &\ the decrease of the 

function x{t, b) is more steep along the t-axis, and changes its monotonic decrease with 
a decrease of value of t at any fixed b changes to an increase at some finite t m (b) > 0, 
which renders the minimal value of x a t given b. Obviously, the decrease of x is n °t 
symmetric with respect to the axes t and b even in the pre-critical region t ~ b ~ 0.01. 
The difference in the behavior of x with respect to t and b is better seen in the strong-5 
limit (t < b). 
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Figure 2: The susceptibility shape function b) for q2D systems and variations of 
t = (0.01, 0.1) and b = (0.01, 0.1). 





b 



Figure 4: The susceptibility shape function x(t, 6) for q2D systems: t = (0.001, 0.1) 
and 6 = (0.1, 1). 

4.2 Strong-5 limit 

For relatively large induction B, the leading terms in the sum (|30|) are 



Q/i . t \ In 6 2a c (l + 6) 
S(k z ,t,b) = — + 6 In : 



+ 6 In 



r(i/6) 
r(i/26) 



(40) 



where 6-independent terms and small term of order 0(bt) have been omitted. This 
expression of S(k z , t, b) is valid for any < b < 1 and does not contain e(k z ). Therefore, 
the result (140!) can be obtained by setting t = (t;ok z ) 2 = in Eq. (1301 and by applying 
properties of the gamma function T(z) [24] . 

In this limiting case, F 3D is related with F 2 d by 



3D 



L z 

~z Fid i 
So 



(41) 



as implied by Eq. ( )29l) . and F 2 d is given by 



F 2D = p55 5(6), (42) 

where S{b) is a short notation of the expression ( 1401 of S(k z ,t,b) ~ 5(0,0,6) in the 
strong-!? limit. 
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The sum (HOI) does not exhibit any singularity. We shall briefly discuss the case e <C b <C 
1 in order to reveal and generalize a preceding result for the diamagnetic moment [25] : 
see also Ref. [16J. For small b we obtain from Eq. (HOI that 

S(b) = bln^ + ^. (43) 
(6-independent terms have been once again omitted). 

Now one may obtain a simple expressions for the free energies of q2D and 3D super- 
conductors. The q2D free energy F 2 d will be 



2D 



pSB 



, V2a c 11B 
m — ; h 



Tlkf 



12B 



(44) 



whereas F 3D is given by Eqs. ( HTj) and (H4"l) . The diamagnetic moment M. 2 d and the 
diamagnetic susceptibility xfia wn l be 



M 



2D 



-pS 



In 



V2< 



11B 



(45) 



and 



X i2D) = (46) 
6B 

respectively. For 3D systems, in accord with Eq. (Hlj) . M.^d = (L Z /£ Q )M.2D, and 
X*f = (W£o)%£f , where .M 2 d and are given by Eqs. (SSJ and (HH}. 

These results are shown in Figs. |3] and HJ For q2D systems, the magnetization shape 
function m = Ai/pS is given by m(t,b) = —dS(0,t,b)/db, where S(0,t,b) is equal to 
S(0,t,b) for T^ksjo-c = 1; see Eq. (|30|) . The function m(t,b) is shown in Fig. [3] for 
* = (0.0001, ...,0.001) < b = (0,001,..., 0.01); q2D systems. When b tends to 0.01, 
the variations of m(b ~ 0.01, t) with t G (10~ 5 , 10 -3 ) are relatively small compared to 
those for t ~ 6 ~ 10~ 3 . At given small £, t ~ 10~ 4 in Fig. [31 x slowly increases with the 
decrease of 6 following the linear low (145 jl . and tends to — In 2/2 « 0.345 for 6 = 10~ 3 
in accord with Eq. ( I45l) ; a result, firstly achieved in Ref. [25J. 

The susceptibility shape function x(t, b) of q2D systems is shown in Fig. HJ for t = 
(0.001, . . . , 0.1) and b = (0.1, . . . , 1). As seen from Fig. HJ the shape function x(t, b) 
remains finite, provided b ^> t even when b tends to zero. Besides, as we have shown 
analytically for b -C 1, in the large- 1? limit this function virtually does not depend 
on t for any fixed < b < 1. In accord with our analytical result Eq. (I46j) . valid for 
t < i) < 1, FigHJ shows that at fixed b, the function x(t, b) is almost constant for 
variations of t under the condition t -C b. At fixed t -C b, however, the variations of 
the function \ip) are substantial, in particular, for b ^ 1. When t tends to zero, the 
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function x(t ~ 0,6) is bounded from below at —11/6, as seen from both Eq. (H6"j) and 
Figll 

4.3 Discussion of the results: application to p-wscve ferromag- 
netic superconductors 

The results for the free energy F, the diamagnetic moment A4, and the diamagnetic 
susceptibility Xdia are very similar to the respective known results for conventional 
non- magnetic superconductors [16] . In particular we point out the dependence of these 
physical quantities on the parameters t and b, describing the departure of thermody- 
namic states from the phase transition line T C (P). For non-magnetic superconductors 
(M = 0), B = H, T C (M = 0) = T c0 is the usual superconducting critical tempera- 
ture at zero external magnetic field, and we reveal the known results for 3D and q2D 
standard superconductors, summarized in the review |16j . 

In p-wave ferromagnetic superconductors, the magnetization M in zero external mag- 
netic field H is different from zero along the whole phase transition line T C (P) and the 
shape of critical temperature T C (P) in zero external magnetic field H strongly depends 
on the magnetization M, as given by Eq. ( )25l) . Thus, except for a very narrow domain 
of the T — P phase diagram above the critical pressure P c , the magnetization above the 
line T C (P) is always large and, hence, for this case, we should consider large values of 
the induction B even when the external magnetic field H is small or equal to zero. The 
important quantity in our consideration is the induction B, because the latter enters in 
the magnetic frequency ojb and the magnetic length ag. Now the magnetic induction 
B plays a role similar to that of external magnetic field H in the theory of diamagnetic 
moment and diamagnetic susceptibility in usual superconductors [El [19]. Therefore, 
for p-wave ferromagnetic superconductors with phase diagrams of the types shown in 
FigJH we should use only those of our results, which correspond to the large- B limit. 
Our results show that both diamagnetic moment and diamagnetic susceptibility do not 
exhibit any singularity and remain finite up to T = T C (P) along the most part of the 
FM-FS phase transition line T C (P). This means that the diamagnetic singularities are 
dumped by the ferromagnetic order. 

The results in the weak-S limit could be valid in a close vicinity of the critical pres- 
sure P c , where the magnetization M of the ferromagnetic phase is very small and the 
criterion for weak-S limit is fulfilled: (H + 4ttM) <C Bq. Then the results in this 
limit will be valid for enough small external field H and M ~ M c = M(T C ) <C Bq. 
In this CclSC, clS mentioned in Sec. IV. A. 1, the diamagnetic moment Ai < will exist 
even for H = and will be proportional to the ferromagnetic moment M ~ M c ; see 
Eq. f l35]) . According to Eq. fl35|) . the 3D superconductor will have a negative total 
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magnetic momentum M tot = A4 + VM at H = provided 



WT > * 1/2 ' (47a) 

According to Eq. (1381) . the total magnetic moment in 2D superconductors will be 
negative, M tot = M + SM < 0, provided 

^ > t (47b) 

When the criteria (I47al) and (I47bj) are satisfied the diamagnetism prevails and the 
overall magnetization of the system is negative. However, under certain conditions, 
some relevant fluctuation contribution of the magnetization vector M may occur, and 
this is an issue which need a study beyond the Gaussian approximation. A reliable 
application of our results in the weak-5 limit could be performed in ferromagnetic 
superconductors, where a line of phase transition of type N-FS exists, i.e., when the 
lines of the N-FM and FM-FS phase transitions of first order meet at some finite tem- 
perature critical-end point (T, P' c ~ P c ) which is connected with the zero temperature 
point (0, P c ) by a second-order (N-FS) phase transition line (Sec. I). 



5 Conclusion 

Introducing an advanced theoretical approach, we have been able to investigate the 
basic properties of the superconducting fluctuations in p-wave ferromagnetic super- 
conductors in zero external magnetic field. For the presence of a strong magnetiza- 
tion due to the ferromagnetic state, the p-wave ferromagnetic superconductors with 
Tjp(P) ^> T C (P) exhibit a diamagnetic behavior, which is typical for usual supercon- 
ductors in the strong- H limit. For this type of ferromagnetic superconductors we have 
demonstrated a form of universality. It is known [16] that at a strong external field 
H, the diamagnetic quantities do not exhibit singularities. Here the same quantities 
undergo the same dumping, i.e. lack of singularities in the strong-I? limit, including 
their values at external field equal to zero. 

In the weak-5 limit the diamagnetic moment and the diamagnetic susceptibility ex- 
hibit scaling singularities with respect to the parameter t of type known from the 
theory of non- magnetic superconductors [I7J. As demonstrated in Ref. [17] for 3D 
geometry, and here for both 2D and 3D superconductors, the scaling amplitudes for 
p-wave ferromagnetic superconductors are quite different from the known scaling am- 
plitudes for nonmagnetic superconductors [16J. The difference is due to the existence 
of ferromagnetic moment M above the FM-FS phase transition line Tps{P)', see Fig. [U 
Therefore, our new results in the weak- 13 limit may have application to p-wave super- 
conductors with T — P diagrams of shape shown in Fig. (TJ, where the FM-FS phase 
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transition line lies below the N-FM phase transition line. Besides, in order to apply 
the weak-S limit results, the ferromagnetic moment M should be enough small. Thus, 
the results for weak-B could be applied only in a close vicinity of the critical pressure 
P c , where the N-FM and FM-FS phase transition lines are very close to each other and 
the ferromagnetic states between them possess a small magnetization M. If the N-FM 
phase transition in this domain of T — P diagram is of first order, as indicated by the 
experimental data, the ferromagnetic fluctuations are suppressed and could not affect 
on the fluctuation diamagnetism. 

The weak-5 limit may be applied to p-wave superconductors containing a N-FS line of 
phase transition. Then the fluctuation diamagnetism in the iV-phase will be described 
by the known formulae [16]. For the lack of ferromagnetic moment in the N-phase 
(M — 0, B — H), in such cases we should consider "weak- if limit [16J. The physics of 
p-wave ferromagnetic superconductors is not limited to the ferromagnetic compounds 
enumerated in this paper and those discovered until now. In future, new substances 
exhibiting p-wave ferromagnetic superconductivity with different shape of the T — P 
phase diagram may be discovered. The theory predicts a variety of possible T — P 
phase diagrams, including both cases with Tp > T s and T s > Tp. In the last case, 
stable pure (non-magnetic) phases are possible [7j. This means that the results for the 
weak- if may have a wider application. 

Except for the location of the phase transition line T C (P) at H = 0, given by Eq. (125]) . all 
results for the diamagnetic quantities p-wave superconductors in Gaussian approxima- 
tion can be obtained from the known results for usual (non-magnetic) superconductors 
by the substitution H — > B. This is a form of universality, deduced in the present 
paper. 

We have used the Gaussian approximation, which is not valid in the critical region [18] 
of anomalous fluctuations. As the critical region of real ferromagnetic superconductors 
with spin-triplet electron pairing is often very narrow and, hence, virtually of no inter- 
est, the present results can be reliably used in interpretation of experimental data for 
real itinerant ferromagnets, which exhibit low-temperature p-wave superconductivity 
triggered by the ferromagnetic order. 
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